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^ ! ABSTRACT 

' Context. First hydrostatic cores are predicted by theories of star formation, but their existence has never been demonstrated con- 

\ vincingly by (sub)millimeter observations. Furthermore, the multiplicity at the early phases of the star formation process is poorly 

constrained. 

Aims. The purpose of this paper is twofold. First, we seek to provide predictions of ALMA dust continuum emission maps from early 
Class objects. Second, we show to what extent ALMA will be able to probe the fragmentation scale in these objects. 
\ Methods. Following our previous paper (Commergon et al. 2012, hereafter paper I), we post-process three state-of-the-art radiation- 

magneto-hydrodynamic 3D adaptive mesh refinement calculations to compute the emanating dust emission maps. We then produce 

■ synthetic ALMA observations of the dust thermal continuum from first hydrostatic cores. 
$^14 Results. We present the first synthetic ALMA observations of dust continuum emission from first hydrostatic cores. We analyze the 

results given by the diff'erent bands and configurations and we discuss for which combinations of the two the first hydrostatic cores 
would most likely be observed. We also show that observing dust continuum emission with ALMA will help in identifying the physi- 
cal processes occurring within collapsing dense cores. If the magnetic field is playing a role, the emission pattern will show evidence 
^ ' of a pseudo-disk and even of a magnetically driven outflow, which pure hydrodynamical calculations cannot reproduce. 

' Conclusions. The capabilities of ALMA will enable us to make significant progress towards understanding fragmentation at the early 

Class stage and discovering first hydrostatic cores. 
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1, Introduction the collapse) issues appear to b e critical to better constrain the 

^— ^ , . , . . star formation mech anism (e.g. jLi et aL I I20HIJ00S et"aD[2QT2l: 

^-^ . It is established that most stars form in multiple systems ISeifried eTaril2012l) 

■ (iDuquen noy & Mayor 1991; Janson et al. 201Z). This indicates 

^ ; a fragmentation process during star formati on, which can be ex- The tremendous combined developments of observational 

' plained by several mech anisms (e.g.. iBodenheimer et al.ll2000l: and supercomputing capabilities allow to study astrophysical 

iMcKee & Ostrikejl2007h. The first picture is to consider the in- processes on scales which were until today unresolved. In par- 

terplay between turbulence and gravity within molecular clouds, ticular, great advances in understanding the star formation pro- 

which can lead to an initial fragm entation prior to the col- cess have been achieved during the past ten years. On the one 

lapse (iHennebelle & Chabriel^ |2008|). In this picture, the stel- hand, thanks to various (sub)millimeter interferometric facili- 

lar Initial Mass Function (IMF) is mainly determined at the ties (e.g., the IRAM Plateau de Bure Interferometer, PdBI, and 

dense core formation stage, the latter under going collapse with- the Submillimeter Array, SMA) and to the Spitzer and Herschel 

out fragmenting into individual objects (e.g., Price & Bate''2007|; space telescopes, much progress has been achieved towards un- 

[Hennebelle & Teyssier 2001; Commercon et al. 2010). On the derstanding the formation and structure of prestellar dense cores, 

other hand, fragmentation may also occur during the collapse and c onstraining the evolu tionary stages of star-forming regions 

of molecular clouds (e.g [Bate & Bonnej |2005|; |Bate| |2012|) (e.g.. iKennicutt & Evansll2012h . On the other hand, numerical 

or within disks that are forn ied because of the conservation models of star formation integrate more and more physical pro- 

of angular momentum (e.g., IWhitworth & StamatellosI 120061: cesses. Among the most important ones, magnetic fields and 

|Commer9on et al.| |2008|). The fragmentation process thus re- radiative transfer appear to shape the collapse and fragmenta- 

mains a matter of intense debate, and in particular, the disk for- tion of pre stellar dense cores (e.g.. IHennebelle & Tevssi^l2008l: 

mation and early fragmentation (i.e., during the early phase of 'Bate'2009'), while their combined feedback dramatically inhibits 

fragmentation in low- and high -mass collapsing dense cores 

Send offprint requests to: B. Commergon (ICommercon et al.l 1201 Ol l2011bh . Unfortunately, there is cur- 

benoit.coinmercon@lra. ens.fr rently no direct evidence of these mechanisms, since observa- 
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tions are not yet able to probe the fragmentation scale in nearby 
star-forming regions. While wide ^ 5000 AU multiple systems 
are often detected around the youngest (C lass 0) protostars (e.g., 
IChen et al ] |2008l:lLaunhardt et al.ll201Qh , the highest resolution 
observations so far with synthesized half power beamwidths of 
0.3" - 1" (e.g., using the IRAM PdBI or the SMA) show a 
lack of close < 2000 AU multiple systems ('M aurv et al]l2010l) . 
However, at a distance of the nearest star-forming regions (e.g., 
140 pc for Taurus), this angular resolution probes only linear 
scales larger than 40-50 AU. In order to probe smaller scales, 
Atacama Large Millimeter/Submillimeter Array (ALMA) obser- 
vations are definitely needed. 

First hyd rostatic cores (FHSC), i.e., first Larson cores 
(iLarsonll 19691) are the first protostellar objects formed during the 
star formation process with typical sizes of a few AU. Al t hough 
their existence is pr e dicted by theory (e . g.. iLarsonl 1 19691: 
Masunaga etaP 1 19981 : iTomida et aP l2010bl: ICommercon et al.l 



2011ah . there still is no strong observational evidence for such 
objects, because FHSCs are deeply embedded within collaps- 
ing cores and their lif etimes are relatively short (at most a 
few th ousand years, e.g.. iTomida et al .1120 1 Oal: [Commercon et al .1 
I2OI2I hereafter Paper I) compare d to the Class phase duration 
(0.1 - 0.2 Myr,"Evansetal."2009^. Several candidate FHS Cs are 
known (e .g., Belloche et al. 2006; Chen et al. 2010; Pined aetaD 
I2OIII: [Chen et all 120121) . but none has yet been confirmed. We 
showed in Paper I that spectral energy distributions (SEDs) of 
star-forming clumps can help in identifying FHSC candidates, 
but are not able to assess the physical conditions within the cores 
and in particular the level of fragmentation. For that, high resolu- 
tion interferometric imaging is needed. It is also unclear to what 
extent FHSCs can be distingu ished from more evolved very-low 
luminosity objects (VeLLOs: ldi Francesco et a n i2007h and sec- 
ond hydrostatic cores (i.e., the prostostars). 

To address this necessity of theoretical pr edictions for the 
appear a nce of early pha s es of sta r formation, iKrumholz et al.l 
(l2007h. 



ISenienov et al.] (l2008h . ICossins et all (120101) . and 

lOff'ner et al.1 ( l2012h presented synthetic ALMA observations of 
dust continuum or line emission. In this paper, we present the 
first predictive dust emission maps of embedded FHSCs as they 
should be observable with ALMA. This study focusing on dust 
continuum is considered as a first step towards FHSC character- 
ization in combination with the results of Paper I. 

The paper is organized as follows. Section 2 presents the 
physical models and the method we use to derive synthetic 
ALMA dust emission maps. Section 3 reports on the results we 
obtain to select the best ALMA configuration and the best re- 
ceiver band in order to observe FHSCs. We discuss the limita- 
tion of our work in Sec. 4. Section 5 presents our conclusions 
and perspectives concerning the future work needed to confirm 
the results obtained as a first step with dust continuum observa- 
tions. 



2. Method 

In this study, we restrict our work to the early stages of the 
star formation process, i.e., the first collapse and FHSC forma- 
tion. As mentioned in the introduction, a lot of progress has been 
achieved in theory and observations to characterize what we 
would expect at those early stages. In the following, we combine 
state-of-the-art tools to produce synthetic ALMA dust emission 
observations. 



2.1. The physical models 

We performed 3D full radiation-magneto-hydrodynamic 
(RMHD) calculations using the adaptive mesh refinement code 
RAMSES ( Teyssier 2002), \ yhich integrates the eq uations of ideal 
magneto-hydrodynamics (iFromang et al.l l2006l) and uses the 
gre y flux-limited-diff^usion ap proximation for the radiative trans- 
fer (ICommercon et al.|[201 id) . We used the same RMHD calcu- 
lations presented in Paper I, which consisted in letting rotating 
(in solid-body rotation) 1 M© dense cores collapse, with initially 
uniform temperature, density and magnetic field. The ratio of 
the initial thermal and rotational energies to the gravitational en- 
ergy are respectively a = 0.35 and jS = 0.045. We ran three 
diff'erent models with the same initial conditions except for the 
initial magnetization, which is parametrized by the mass-to-flux 
to critical mass-to-flux ratio ji = (Mo/0)/(Mo/0)c. The three 
models are depicted as follows : MU2 model (strong magnetic 
field, 11 = 2), MUlO model (intermediate magnetic field, yu = 10) 
and MU200 model (quasi-hydro case, jd = 200). These three 
models are representative of the diversity of FHSCs and of their 
environments (disk, pseudo-disk, and outflow) that are predicted 
by the theory. The difl'erent physical structures and FHSC life- 
times found in the three models are summarized in Table [T] 
In the MU2 model, only one FHSC is formed, surrounded by 
a pseudo-disk, and an outflow has been launched. The MUlO 
model does not fragment either and results in a system composed 
of a disk, a pseudo-disk and an outflow. The MU200 model 
has classical features of hydrodynamical models, in which rel- 
atively large disks (~ 150 AU) are formed and subsequently 
fragment. Readers are referred to Paper I for a thorough de- 
scription of the difl'erent models and their limitations. Note that 
the MU2 and MUlO models are more representative of the ob- 
served magnetization l evel in star forming re gions (i.e., ^ 2-3, 
iFalgarone et al.ll2008l: [Crutcher et al.ll2010h . even though deter- 
mining magnetic field strength remains challenging. 

In the following, the calculations are post-processed at the 
same times as in Fig. 2 of Paper I. These are representative of 
the characteristics of the three models, i.e., + 0-78 kyr for the 
MU2 model, + 1-81 kyr for the MUlO model, and to + 3.26 
kyr for the MU200 model (where corresponds to the FHSCs 
formation in each model, i.e. 50.4 kyr for MU2, 35.7 kyr for 
MUlO , and 35.4 kyr for MU200). 



2.2. Dust emission models with RAMDC-3D 

We used the same interface as that presented in Paper I, 
which couples the outputs of the RMHD calculations, done 
within RAMSES, to the 3D radiative transfer code RADMC-3eQ. 
In this interface, we assumed that gas and dust are thermally 
coupled (rdust = T'gas), which is a valid ap proximation give n 
the high density within the dense cores (e.g. iGalli et al.|[2002l) . 
We also assumed that the gas tem perature T^as computed in 
the RMHD calculations is correct (ICommercon etal 201 Id: 



Vavtet et al.ll2Q12i). W e used the low temperature opacities of 
Semenov et al. for a model in which dust is made of ho- 

mogeneous spheres with a "normal" (Fe/Fe -h Mg = 0.3) silicate 
composition. 

The dust thermal continuum emission maps were computed 
on a square box with physical size Ax = 1300 AU and a res- 
olution 6x = 2.54 AU. Since models were assumed to be at 
a distance Z) = 150 pc, this translates into an angular extent 
A'& = 8.67" and an angular resolution Si9^ = 0.017". Four diff'er- 



www . ita . uni-heidelberg . de/~dullemond/sof tware/radmc- 3d/ 
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Table 1. Summary of the physical structures and FHSC lifetimes found in the three models. The last four columns indicate whether 
these structures are observed (cross) or not (dash) in the synthetic ALMA observations presented in Fig. [3] to [6l 



Model 


Physical structures 


FHSC lifetime 


Structures in synthetic observations 
61 = 0° 6* = 45° 6* = 60° 6* = 90° 


MU2 


one object 
pseudo-disk 
outflow 


1.2 kyr 


X 
X 


X 
X 


X 
X 


X 
X 


MUlO 


one object 
pseudo-disk 
compact disk 
outflow 


1.8 kyr 


X 
X 
X 


X 
X 
X 


X 
X 
X 


X 
X 

X 


MU200 


five fragments 
disk 


> 3 kyr 


X 
X 


X 
X 


X 
X 


X 



Table 2. Characteristics of the ALMA bands used, denoted B. 
Listed are the central frequency vq, full bandwith Avq (not to 
be confused with the bandwidth Ay = 8 GHz accessible for a 
given observation), the field-of-view at the band center, and the 
availability of the band at the Early Science stage. 



B 


vo [GHz] 


Avo [GHz] 


FoV ["] 


Early Science 


3 


100 


32 


63 


Yes 


4 


144 


38 


44 


No 


6 


243 


64 


26 


Yes 


7 


324 


98 


19 


Yes 


8 


442.5 


115 


14 


No 


9 


661 


118 


9.5 


Yes 



Table 3. Characteristics of the full ALMA configurations used, 
denoted C. Listed are the minimum and maximum baselines, and 
synthesized beam major and minor axes for bands 3 and 9. 



c 


^min [m] 


^maj [m] 


O3 ["] 


O9 ["] 


5 


15 


390 


2.47x2.12 


0.37 X 0.32 


10 


20 


924 


1.07x0.99 


0.16x0.15 


15 


25 


1814 


0.48 X 0.43 


0.07 X 0.06 


20 


49 


3699 


0.24 X 0.22 


0.04 X 0.03 



ent viewing angles are used in the following : ^ = 0° (equatorial 
plane seen face-on), 6 = 45°, 6 = 60°, and = 90° (edge- 
on view, perpendicular to the rotational axis). These model dust 
emission maps were computed for six ALMA bands (3, 4, 6, 7, 
8 and 9), including the four that are in use in the Early Science 
stage (see table O, over a total bandwidth Ay = 8 GHz centered 
on the band's central frequency vq. This was done by computing 
a set of ten emission maps, one for each of ten 800 MHz sub- 
bands, and averaging them. The resulting map gives the mean 
specific intensity (ly) of the thermal dust emission over the band, 
in erg cm"^ s"^ Hz"^ sr"\ and is then converted to a brightness 
temperature map in K vi43 



2.3. Producing synthetic ALMA observations 

The GILD AS software packag^E developed and maintained 
by IRAM, is primarily intended for the reduction and analysis of 

^ The factor 10~^ comes from the CGS to SI conversion of the specific 
intensity, as the ALMA simulator assumes input brightness distributions 
to be in K. 

^ iram.fr/IRAMFR/GILDAS/ 



observational data acquired via the IRAM instruments (single- 
dish 30-m radiotelescope at Pico Veleta and 6-antenna array at 
IRAM PdBI). It also includes a full ALMA simulator, with up- 
to-date array configuration^, which was primarily developed to 
assess the impact of the ALMA compact array (AC A) on the 
imaging capabilities of ALMA ( Petv et al. 2001; Tsutsi mfiet al.l 
12004]) . 

To produce synthetic ALMA observations from our dust 
emission models, maps are converted to the GILDAS data format 
(GDE) and projected on a fixed sky position (ao, ^0) = (0°, -23°) 
such that sources transit at the zenith of the array. The field-of- 
view in each band (see tableO, given by 1 .22Ao/d with J = 12 m 
the antenna diameter and Aq the central wavelength of the band, 
shows that a single pointing sufl^ices to map the emission in all 
cases. 

ALMA consists of two sub-arrays, the fifty 12-m antennas 
and the twelve 7-m antennas of the ACA. The array operation 
will consist in constantly moving antennas around, so that no 
two observations may be obtained with exactly the same con- 
figuration. We do not consider the compact array and focus on 
imaging the cores' thermal emission with ALMA only. To keep 
the number of simulations down to a manageable number, we 
select four "typical" configurations (out of the 28 representative 
configurations implemented in GILDAS), whose properties are 
listed in table [3] Since we make one simulation per model, per 
configuration, per inclination angle, and per band, the total num- 
ber of simulations is thus 3x4x4x5 = 240. In all cases, model 
sources are observed for a total integration time of 18 minutes 
centered on the transit. 

The ALMA simulator is a full pipeline, includi ng the decon - 
volution of dirty images via a CLEAN algorithm (IClarklll980h . 
Outputs can thus be compared directly to input images smoothed 
to the same resolution, which the simulator also provides. 

3. Results 

3.1. A necessary compromise 

In our models, thermal dust emission at millimeter and sub- 
millimeter wavelengths mostly comes from structures that are 
larger than the FHSC (i.e., from the disk and the pseudo-disk). 
Yet, to discriminate between diflTerent magnetization levels fi, 
one needs to resolve the fragmentation scale of a few AU. This 
latter constraint means that only extended configurations of the 
array will be able to provide a suflScient angular resolution, in 
the few tens of milliarcseconds range. However, this comes with 

^ We used release apr 1 Ih of GILDAS. The latest releases also contain 
Early Science Cycle 1 configurations of the array. 
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Fig. 1. Left : Probability density function from the distribution of samples in the u - v plane for the 4 configurations used here 
(numbers 5, 10, 15 and 20, indicated next to each curve) as a function oiu-v radius J^v = Vw^ + v^. The relevant parameters are 
the source position and duration of observation, which are described in section [23l The small vertical lines indicate the minimum 
baseline for each configuration. Also marked, with dashed lines, are characteristic baselines duv = cD/(27iyod) corresponding to a 
physical size J = 10 AU at a distance D = 150 pc, for the six frequencies vq. Right : Power spectra of the input maps at 100 GHz 
(solid lines) and 661 GHz (dashed lines) for ^ = 0°. The wavenumber axis has been rescaled to match the w - v radius of the left 
plot. The vertical black lines mark the positions of the 10 AU scale at 150 pc at these frequencies (solid for 100 GHz, dashed for 
661 GHz, reported from left plot). 



an increase of the minimum baseline length, so that large-scale 
emission is lost due to the central hole in the visibility (Fourier) 
space, also called u-v plane. 

To make this idea more quantitative, ¥ig.\l\(left) shows the 
distribution of visibility samples as a function of radius in the 
u-v plane for the four configurations. Marked on this graph 
are the baselines duv = cD/ilnvod) corresponding to a physical 
size J = 10 AU (roughly equal to the size of fragments in the 
MU200 models) at a distance Z) = 150 pc, for the six central 
frequencies vq. On the right plot of Fig. [T] we display the power 
spectra of the face-on input maps at 100 GHz (band 3) and 661 
GHz (band 9), with the wavenumber axis k rescaled to match 
that of the u - V radius on the left plot. This is done by noticing 
that the largest wavenumber in the power spectrum corresponds 
to the pixel physical size Sx. 

The fragmentation in the MU200 model appears as the ex- 
cess power at intermediate scales compared to MU2 and MUlO. 
This excess is clearly seen for 300 m < duv ^ 2000 m at 100 
GHz, with duv = 2000 m roughly corresponding to the 10 AU 
scale at 150 pc. As far as the four configurations used here are 
concerned, this range of baselines is beyond C = 5 and best 
probed by configurations C = 15 and C = 20. However, both of 
these configurations have a larger central hole than C = 5 and 
C = 10, and therefore should lose a larger amount of flux. This 
flux loss becomes larger at higher frequencies, since duv ^ 1 /^o 
implies a global compression of the sources' power spectra to- 
wards smaller baselines. 

For a global view on the emission loss at large scales, Fig.[2l 
shows the ratio between the observed flux to the model flux as 
a function of frequency and array configuration. As expected, 
the most compact configurations (C = 5 and C = 10) allow to 
recover essentially all of the flux in bands 3 to 7. Only in band 
8 and 9 we see a ~ 50% drop in the received flux for the non- 
fragmenting models MU2 and MUlO. Figure [2l also confirms 
that flux loss increases with frequency in configurations 15 and 
20, most (~ 80%) of the flux being lost already at 144 GHz in 
configuration 20. Overall, it appears that observing in bands 3 



and 4 with configuration 15 may provide the best compromise, 
as the flux loss is very limited (less than 8%) and according to the 
left plot of Fig.[Tl it should be able to resolve the fragmentation 
scale, at least in band 4. 

3.2. Emission maps 

We show, in Fig.[3]to[6l the brightness distribution maps out- 
put by the ALMA simulator in bands 3 and 4, for configurations 
15 and 20. In each figure, four inclination angles are shown, 
i.e., 6 = 0°,45°,60°, and 90°. Also shown are the sensitivity 
limits 3crs in these two bands, given by the ALMA Sensitivity 
Calculatofl The most extended configuration and the largest fre- 
quency best probe the fragmentation scales (C = 20, B = 4 in 
Fig.O. However, as we discussed above, a significant amount of 
the flux emanating from the extended emission is lost in config- 
uration 20. On the other hand, the fragmentation scale is barely 
resolved with C = 15 and B = 3. A better compromise between 
high resolution imaging and flux recovering is found using ei- 
ther C = 15 and 5 = 4 or C = 20 and 5 = 3. Overall, there is a 
clear distinction between the two magnetized models (MU2 and 
MUlO) and the quasi-hydro MU200 model. As expected theo- 
retically, as soon as a magnetic field is taken into account, with 
field strengths i n the range of what is act ually observed (e.g ., 
iFalgarone eta T 2008; ' Crutcher etaT]l2010l : iMaurv et al.ll2012h . 
the picture changes dramatically. 

Interestingly, many features that are directly probing the 
physical conditions can be observed with the dust emission al- 
ready (see Table [B. In the MU200 model, the fragmentation is 
resolved for inclination angles 6 < 60°, and only the compact 
emission of the disk is observed in the edge-on view. In con- 
trast to that, the dust emission is much more extended in the 
magnetized MU2 and MUlO models, and corresponds to den- 
sity features that are typical of magnetized collapse: the pseudo- 
disk and the outflow. In both cases, the most powerful emission 

^ The weather conditions used are that of the first octile, i.e., 0.47mm 
of precipitable water vapor above the instrument. 
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Fig. 2. Ratio i^obs/^modei of the observed flux to the model flux as a function of frequency, for the three difl'erent models (MU2, 
MUlO, and MU200), four array configurations (C = 5, 10, 15, 20), and four inclinations (0 = 0°, 45°, 60°, 90°). Magnetization level 
jd increases from left to right, and inclination angle 6 increases from top to bottom. 



comes from the central FHSC (point-like source,yellow area). In 
the MUlO model, the disk emission is observed (blue regions) 
for inclination angles 6 < 60°, but it remains relatively small in 
comparison to the pseudo-disk one. The pseudo-disk emission 
features are very similar in the MU2 and MUlO models for in- 
clination angles 6 < 60°. In the edge-on view (6 = 90°), the 
emission from the pseudo-disk is clearly observed in the MU2 
and MUlO models, and is much more extended than that of the 
disk in the MU200 model. The pattern of the pseudo-disk ob- 
served with C = 20 and B = 4 in the MU2 model may lead to 
confusion since it has the shape of a flared disk. Additional line 
emission observations are then needed to distinguish between a 
disk and a pseudo-disk in that case. Last but not least, the outflow 
is also observed in the MUlO model in Fig. [3] to [5l with an extent 
of ~ 4'' along the polar axis, corresponding to ~ 300 AU. This is 
not surprising since the outflow tends to transport more mass in 
the case of a low er magnetization (Henne belle & Fromang 2008, 
and see Fig. 2 in lCommercon et ani2012l) . 

We investigate thoroughly in appendix (A] the efl'ects of noise 
(thermal and atmospheric phase) and pointing errors on the syn- 
thetic observations. The main conclusions draw from the noise- 
free case remain unchanged, except that more flux is lost in noisy 
observations. We refer interested readers to appendix lAl 



4. Discussion 

One limitation of this work is that w e used a unique dust 
opacity model from lSemenov et al" ] (l2003h . in which dust grains 
can grow, so that we do not predict the variation of the dust 
emission with difl'erent dust grain properties. The dust properties 



(size, composition, and morphology) within dense cores are still 
very uncertain, but there is nevertheless theoretical and observa- 
tional evidence of dust grain growth with i n dense cores (e.g. , 
lOrmel et al.l l2009l: ISteinacker et"aDl2010l) . lOrmel et al.l (l201lh 
compared the dust graiii opacit ies they go t from the various grain 
models of lOrmel et aP (l2009l) with the lOssenkopf & HenningI 
(19^ opacities model and found that the opacity, in the sub- 
millimeter wavelength range, only varies by a factor of a few 
between all the difl'e rent models. We also checked that the 
ISemenov et al.l (l2003l) opacities w e used are similar to the from 
the'Ossenkop f & HenningI (1 19941) ones by factors close 1 , so that 
the choice in the opacity will not change the picture. In addition, 
we show in Paper I (Fig. 4) that the radius at which the optical 
depth equals unity does not depend strongly on the opacity for 
(sub)millimeter wavelengths. We thus speculate that our predic- 
tions are relatively robust given the small variations in optical 
depth and in opacity in the (sub)millimeter wavelength range. 

The second main limitation comes from our idealized initial 
and boundary conditions, which do not account for the dense 
core environment (e.g., eventual mass accretion on the core 
while it collapses) and initial turbulence, whereas the latter can 
potentially modify the magnetic braking and thus the fragmen- 
tation properties during the collapse (e.g., for higher mass dense 
cores, ISeifried et al]l2012h . Turbulence, howeve r, is observed to 
be siib- to trans-sonic in low mass dense cores (iGoodman et al.l 
119981: lAndre et al.ll2007b . so that it will not change dramatically 
the outcome of the magnetized dense core collapse. In addition, 
magnetic fields dominate the dynamics and inhibit the fragmen- 
tation even more when combined with radiative transfer because 
of the energy released from the accretion shock at the FHSC bor- 
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Fig. 3. Brightness distribution maps, in Jy/beam, output by the ALMA simulator in band 5 = 3 (100 GHz) and array configuration 
C = 15. Magnetization level ju increases from left to right, and inclination angle increases from top to bottom. Contours show the 
3crs sensitivity limit in this band, as given by the ALMA Sensitivity Calculator. Here os = 14.55 yuJy. The synthesized beam is 
shown in the bottom left corner of each plot. For a given 6, color scales are identical across all columns. 
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Fig. 4. Same as Fig.[3]but for 5 = 4 and C = 15. Here, crs = 16.05 yuJy. 
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Fig. 5. Same as Fig.[3]but for 5 = 3 and C = 20. Here, crs = 14.55 yuJy. 
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Fig. 6. Same as Fig.[3]but for 5 = 4 and C = 20. Here, as = 16.05 yuJy. 
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der (ICommercon et al.ll2010ll2011bh . We thus conclude that the 
three models presented here are representative of the variety of 
FHSCs that can be formed with various initial conditions. 

Finally, all the analysis and the ALMA synthetic map calcu- 
lations have been done for objects placed at a distance of 150 
pc. The best compromise between flux loss and angular reso- 
lution may thus be altered for objects at significantly diff'erent 
distances, but the method remains unchanged. 



5. Summary and perspectives 

We present synthetic dust emission maps of FHSCs as they 
will be observed by the ALMA interferometer, using state-of- 
the-art radiation-magneto-hydrodynamic models of collapsing 
dense cores with diff'erent levels of magnetic intensity. We post- 
process the RMHD calculations performed with the RAMSES 
code using the 3D radiative transfer code RADMC- 3D to produce 
dust emission maps. The synthetic observations are then com- 
puted using the ALMA simulator within the GILD AS software 
package. 

We show that ALMA will shed light on the fragmentation 
process during star formation and will help in discriminating 
not only between the diff'erent physical conditions, but also be- 
tween the diff'erent models of star formation and fragmentation. 
The forthcoming facility will yield interferometric observations 
of FHSC candidates (selected using SED data, see Paper I) in 
nearby star-forming regions with sufficient angular resolution to 
probe the fragmentation scale. We also show that the intensity 
of the magnetic field in this early phase of protostellar collapse 
will also be assessed, since we do see clear morphological differ- 
ences between the magnetized and non-magnetized models, i.e., 
pseudo-disk and outffow in the former case versus disk and frag- 
mentation features in the latter case. We stress that due to the un- 
precedented sensitivity and coverage of the instrument, this will 
be achievable in a very short time (18 minutes), so that many 
FHSC candidates may be observed in a single observing run. 

Reaching good enough angular resolution to probe the frag- 
mentation scales and limiting the loss of large-scale emission re- 
quires a compromise, which is achieved using frequency bands 
3 and 4 and the relatively extended configurations 15 and 20. 
We also investigate the effect of noise on the interferometric ob- 
servations and show that in typical conditions (see Fig. lA.lK 
ALMA will still be able to reveal fragmentation. We note that 
the effect of atmospheric phase noise can be efficiently reduced 
using water vapor radiometers on the ALMA antennas. 

We do not discuss the impact of the ALMA Compact Array 
(AC A) on our simulated observations. Its purpose is to pro- 
vide measurements for the large-scale emi ssion and improv e the 
wide-field imaging capabilities of ALMA (iPetv et al.ll200lh . but 
it is not meant to resolve fragmented molecular cores. In our 
models, this requires a 0.5" angular resolution (see e.g., the case 
B = 3 and C = 15 on Figs. [T] and [3]), but the longest baselines 
accessible to AC A are ~ 100 m, so that the best angular resolu- 
tion available, at the high-frequency end of band 10 (950 GHz), 
is 0.8". However, AC A could help in recovering some of the lost 
flux from the disk, pseudo-disk and outflow. 

Our work is currently limited to the dust continuum emis- 
sion, which cannot yet provide robust means to discriminate be- 
tween FHSC and second hydrostatic cores (SHSC) and to con- 
clude on the nature of VeLLOs. Further work including molec- 
ular line emission calculations is thus warranted to better probe 
the physical conditions (density, temperature, etc..) in observed 
collapsing cores, for instance to disentangle between the disk 



and the pseudo-disk, which should harbour different line pro- 
files (rotation dominated versus infall dominated). Line emission 
predictions are thus the obvious next step towards FHSCs char- 
acterization. 
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Appendix A: Noisy observations 

The simulations presented in the main body of the paper 
were performed in the unrealistic case of noiseless observations. 
To truly assess ALMA's ability to resolve fragmented cores, an 
extension of this study to noisy observations is required, taking 
into account the different causes of noise : pointing errors, ther- 
mal noise and phase noise. 

We start with an error-free simulation of the MU200 model 
viewed face-on (6 = 0°), in band B = 4, with configuration 
C = 20. All parameters are identical to those of the general set 
of simulations, but we add zero- spacing with a single-dish ob- 
servation, so that in this case all of the flux is recovered (0.1090 
Jy) instead of 96% of it (0. 1046 Jy) with ALMA only. 

The following errors, whose results are summarized in 
Table lA.ll may then be added, separately or simultaneously, to 
the simulated observations : 

► 0.6" random pointing errors cTq,, which correspond to ~ 15% 
of the primary beam width at 144 GHz. They may be applied 
to the ALMA antennas alone or to both ALMA and single-dish 
measurements. 



Table A.l. Efl'ect of the difl'erent noise types. The error- free sim- 
ulation is that of a combined B = 4 observation of the face- 
on MU200 model with ALMA in configuration C = 20 and 
with a single-dish. Pointing errors (cTq,) and thermal noise (ctt) 
may be applied to neither instrument ("No"), to the interferom- 
eter only ("ALMA"), or to both interferometer and single-dish 
("Both"). Phase noise is specified via the rms atmospheric phase 
on a 300-m baseline, which can take the values 30° or 45°, and 
can be corrected via water- vapor radiometers ("WVR"). Shown 
are the fluxes S of the output maps and the median fidelities 
Fi, F3 and Fio, on pixels whose intensities in the model image 
are higher than 0.3%, 1%, 3% and 10% of the peak, respectively 
(see section lAt. 
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► Thermal noise ctt. In band 4, the ALMA Sensitivity 
Calculator suggests using the 6* octile for the atmospheric con- 
ditions, which corresponds to 2.75 mm of precipitable water 
vapor and a zenith opacity r = 0.05 7. The band 4 r eceive r 
noise is set to 7144 = 40 K following lAsayama et al.l (l2008l) . 
Regarding the single-dish measurements, we use a system tem- 
perature Tsys = 100 K. The simulator allows for setting none, 
either or both ALMA and single-dish thermal noises. 

► Atmospheric phase noise cr^. The ALMA simulator allows 
for a turbulent atmospheric screen to pass over the interferome- 
ter, distorting the waveplanes and causing phase errors. This 2D 
phase screen is characterized by a second-order structure func- 
tion that is the combination of three power-laws in three spa- 
tial ranges, scaled so that the rms phase diff'erence for a 300-m 
baseline takes a specifi c value. In our ca se, we chose this value 
to be either 30° or 45° (iPetv et al.ll200lh . Situated at an altitude 
z = 1000 m, the screen passes over the instrument at a windspeed 
w = 10 m.s"^ Calibration, which is done every 26 seconds using 
a calibrator 2 degrees away from the source perpendicularly to 
the wind direction, may include the use of water vapor radiome- 
ters (WVR), which directly measure the amount of precipitable 
water vapor along the line of sight in the atmosphere above each 
antenna. 

Table lA.ll gives simulation results associated to the various 
noise situations considered : fluxes in the output (deconvolved) 
maps, and median fidelities on pixels whose intensities in the 
model image are higher than 0.3%, 1%, 3% and 10% of the 
peak. The fidelity is basically the inv erse of the relati ve error 
between the output map and the model (iPety et al .11200 ih . so that 
the higher the fidelity, the better the reconstruction. What is ap- 
parent is that atmospheric phase noise has the strongest impact 
on the reconstruction process, with a third of the flux being lost 
in the worst-case scenario, and fidelities dropping to a few (30%- 
50% relative error on the output maps). The use of WVR is a def- 
inite plus in this situation, as relative errors then drop to a little 
over 10%. However, even in the worst possible situation consid- 
ered here, and without WVR correction, ALMA still is able to 
uncover fragmentation, as Fig. lA. II shows that all five fragments 
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Fig. A.l. Brightness distribution maps, in Jy/beam, for two cases 
considered in section (A] The top plot shows the case of an error- 
free observation at 144 GHz with ALMA in configuration C = 
20, combined with single-dish. The bottom plot shows the case 
of an observation with the same combination of instruments, but 
with 0.6" pointing errors on all antennas, thermal noise as de- 
scribed in the text, and uncorrected atmospheric phase noise with 
a 45° rms phase diff'erence on 300-m baselines. Contours on both 
plots correspond to the 3crs sensitivity limit given by the ALMA 
Sensitivity Calculator, where os = 11 Al yuJy for the typical at- 
mospheric conditions in band B = 4. 



in the model map are well observed above the noise level in these 
typical conditions. 



